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The angular power spectrum of the Cosmic Microwave Background has been measured out to 
sufficiently small angular scale to encompass a few acoustic oscillations. We use a phenomenological 
fit to the angular power spectrum to quantify the statistical significance of these oscillations and 
discuss the cosmological implications of such a finding. 
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During the past few years there has been unprecedented 
progress in the detection and characterization of fluctua- 
tions in the Cosmic Microwave Background on small an- 
gular scales. The TOCO, BOOMERanG and MAXIMA 
experiments presented unambiguous evidence for a 
preferred scale (a "peak" ) in the variance of temperature 
fluctuations as a function of scale, (AT)^ = l{£-\-l)Ct/ A-k 
(where Ci is the angular power spectrum at a scale 
£ ~ ISO^/fl). This is a striking result, having been 
predicted in 1970 from simple assumptions about scale- 
invariance and linear perturbation theory of General Rel- 
ativity. During the last six months the possible presence 
of peaks and troughs in (AT)^ has been reported by the 
BOOMERanG [|, DASI g and MAXIMA [g experi- 
ments. This harmonic set of features in the power spec- 
trum could further confirm that the origin of structure 
was due to a primordial set of fluctuations which set the 
cosmological plasma "ringing" at very early times. It is 
the purpose of this short note to quantify the confldence 
with which one can claim that there are acoustic oscilla- 
tions present in the current data. 

We shall flrst outline the argument for why we expect 
oscillations in the power spectrum of the CMB. Let us 
restrict ourselves to the radiation era, where the photons 
and baryons are tightly coupled. Its suffices to consider 
the density contrast in radiation, 6^ and (in the syn- 
chronous gauge) the trace of the spatial metric pertur- 
bations h = h\. The first order Einstein and continuity 
equations in Fourier space are: 

••1-6 ■■ k"^ 2 •• 

h+ -h + —5^1 = {} and 5^, + —5^ + -h = Q (1) 

where rj is conformal time and k labels the Fourier com- 
ponent. Let us restrict ourselves to adiabatic initial con- 
ditions. In this case one has that 5^ = — One can 
solve Eq. |] on large scales fcr; <C 1 to find two solutions 
6-y = Ak?7 , -Bk?/"^; if these are setup early in the radi- 
ation era, the growing dominates very rapidly and one 
can to as excellent approximation set i?k = 0. On small 
scales one can solve the system using a WKB approxima- 
tion to find S-y (X cos{kr]/V^),sm{kri/\^). Matching the 
large scale solution to the small scale solution one finds 
that 5j = Ai^{kr])'^ cos{k?] / \/3 + (j)) . The acoustic oscil- 
lations in the AT^ will be primarily the power spectrum 
of S-y at recombination, rj^ projected until today: 
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FIG. 1. Data from BOOMERANG, DASI and MAXIMA. 
Overplotted, our best fit from the class of models described 
by Eq. 3. 

AT/ (X (|Ak|2> cos^(^ + U=kxivo-v,) (2) 

where x(^) is the conformal distance corresponding to 
the coordinate distance r. Thus, the oscillations in the 
baryon-photon plasma will lead to a set of peaks and 
troughs in the angular power spectrum. 

How general is this argument? As stated above we are 
considering primordial, passive, perturbations which are 
in effect initial conditions for the evolution of the coupled 
baryon/photon plasma responding to gravity and pres- 
sure 1^. The most general class of such perturbations 
was classified in |^ where, in the context of the current 
menagerie of matter candidates, it is believed that there 
are five degrees of freedom, possibly correlated among 
each other (corresponding to a 5 x 5 symmetric matrix 
of initial conditions). Although we looked at the spe- 
cific case of adiabatic perturbations, the argument fol- 
lows through for all other pure primordial perturbations. 
By pure perturbations we means perturbations in which 
one only picks one of these degrees of freedom to be non- 
zero. The key feature is that the large scale solutions to 
Eq. Ill have two modes, one of which is decaying and very 
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rapidly becomes subdominant. For example for density 
isocurvature perturbations a dilferent phase, 4> will be 
picked out such that the positions of the peaks will be 
out of phase with regards to the adiabatic models. If one 
considers sums of initial condition such as in [|, then it 
is conceivable that the combination of power spectra will 
interfere in such a way as washout the oscillations. 

Clearly, the presence of oscillatory features would be 
strong evidence that the structure was seeded at some 
early time, before recombination. With the current CMB 
data it has been suggested that we are already seeing ev- 
idence for such features. One of the reasons for such a 
claim is that primordial, passive models supply a good 
fit to the measured angular power spectrum and, as ar- 
gued above these models have acoustic oscillations. The 
concern is that, all models which have been compared to 
the data have oscillations in the Ce and one therefore has 
not strictly tested for the presence of these oscillations. 
We propose to do this in the following: construct a phe- 
nomenological fit to the data points that can smoothly 
interpolate between a model with no oscillations to one 
which has oscillations of a well defined frequency and 
phase. An analogous approach was used in |^ to quan- 
tify the significance of the presence of a peak at ^ ~ 200. 
The parametrization we shall use is of the following form: 

ATI ^ Ae'^^^ +Bcos^[a{£-£p) + <j)]+C (3) 

This is a seven parameter family of models and we can 
justify it as follows. The goal is to detect the presence 
of prefered frequency in the AT^ so the parameters we 
are ultimately interested in are a and B. However we 
know that there is a well defined peak in the data with 
a well defined width, and the signal to noise of this peak 
is sufficiently high that it will dominate any statistical 
analysis; i.e. the width of the peak, a will tend to peak 
at a frequency of order tt/ct. Given that we wish to be 
conservative we consider a part of the curve correspond- 
ing to the first peak (characterized by A, ip and a) and 
marginalize over these parameters. Hopefully in this way 
we decrease the statistical weight of the peak. Finally we 
introduce an offset, C which allows us to interpolate be- 
tween a fiat and oscillatory curve. 

Given that we are interested in the behaviour of the 
power spectrum in the regime where acoustic oscillations 
will dominate, i.e. on scales larger than the sound hori- 
zon at last scattering, we shall not include the COBE 
data set. We restrict ourselves to three data sets, 
BOOMERanG @, DASI §] and MAXIMA g. We shaU 
minimize the fitting function of Eq. ^ with regards to the 
three data sets using a standard . A few comments are 
in order. We do not consider more refined parametriza- 
tions of the band-power distribution functions jsj , such as 
the off-set log-normal or the skewed approximation. We 
marginalize over the calibration uncertainties quoted in 
[g-Q . Furthermore the beam uncertainties are taken into 
account by adding them in quadrature to the noise co- 
variance matrices of each experiment. All these approxi- 
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TABLE I. The mean values and 95% errors for the a and 
B as defined in Equation 3 and the corresponding for the 
best fit models (where Nd is the number of parameters). See 
text for description of the different data combinations. 

mations may introduce a modest degree of uncertainty in 
our results but do not change the essential conclusions. 

In table Q we summarize our results for a, B and the 
of the best fit model to each subset of data. The dif- 
ferent combination and subsets of the data we consider 
are: data from all three experiments ( "AU" ) and for all 
three data experiments discarding all points with £ < 400 
("All>4oo"), the data from each individual experiment 
("BOOMERanG", "MAXIMA" and "DASI"), data from 
DASI discarding the point at £ = 553 ("DASI-i"). We 
have also generated a mock realization of the best fit adia- 
batic model to the data with corresponding variance from 
sampling and noise from the combination of the three ex- 
periments ( "Mock" ) and finally the same realization but 
discarding aU points with £ < 400 ("Mock>4oo"). We 
present the mean values and the 95% confidence limits 
from the integrated likelihoods. We find that the data 
does seem to pick out a favorite frequency of oscillation 
in the AT^ of a = 1-1° § ij corresponding to an interpeak 
spacing of S£ — tt / a = 286+44'^. It is interesting to note 
that even removing the points that lie in the region of 
the I''* peak, the detection persists albeit with larger er- 
ror bars. The three experiments detect similar values 
of a with varying confidence regions. We should note 
that the likelihoods in a are extremely skewed and in 
some cases actually have two local maxima. For exam- 
ple the maximum of the likelihood for BOOMERanG is 
a = 0.01 while for MAXIMA there is a local maximum 
a.t a ~ 0.011. Note also the importance of the point at 
£ — 553 in the DASI data. If we remove this point from 
our analysis, the DASI confidence region for a enlarges 
considerably. We have listed the values of the for the 
corresponding best fit models along with the number of 
data points used in each case. All of the x^ are reason- 
able although correlations between the data points may 
lead to a smaller number of degrees of freedom than sim- 
ply Nu — Np where Np = 6 is the number of parameters. 
The Mock data sets lead to similar values of a and cor- 
responding confidence limits. 

Our analysis indicates that there is a marginal pres- 
ence of oscillations in the measured AT^ (at the 2 — a 
level) within the context of the phenomenological model 
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described by Eq. 3. This caveat is important. Although 
we have attempted to justify the functional form of our 
model using rough general arguments, it is conceivable 
that one could construct other models which interpolate 
between oscillatory and non-oscillatory behaviour and 
which reduce or enhance the significance of detection of 
a. For example, if we consider the All>4oo combination 
of data and add a term of the form D£ to Equation 3., 
the mean value of a is still 0.01 but now a — lies within 
the 95% confidence region. We have, of course, done this 
by adding yet another parameter. However this does not 
exclude the possibility that there are models with fewer 
parameters and no oscillations that may provide a better 
fit to the data. 

Let us pursue the implications of the constrain on a 
within the context of models which predict oscillations. 
The above discussion leads us to note that the spacing 
between peaks S£ can be used to set a constraint on the 
angular-diameter distance. As noted in the physi- 
cal peak separation is set by the sound-horizon (which 
we know from atomic physics). Combined with the our 
knowledge of time of recombination and d£ (which is ef- 
fectively half the angular scale subtended by the sound 
horizon at recombination) we can constrain the geometry 
of the universe. Efstathiou and Bond jl^ have proposed 
a convenient parametrization using the "shift" parameter 
R = 2y/QK/^M/xi'>lo - r]*) where fi/f = 1 - f^A - 
and Qm {^a) are the fractional energy densities of mat- 
ter (cosmological constant) today. Indeed one has that 
5£ ~ 290/i? and one can obtain a likelihood for R: one 
finds that R = LOlt^ j?. This can be reexpressed into a 
constrain on 17 = 1 — fix if we marginalize over Q\ to ob- 



discussions and advice. 



tain r2 ~ 1 



Note that this constraint is independent 



from the constraint obtained from the assumption of adi- 
abaticity and the position of the first peak, £p ~ 215/i? 



where one finds R 



1.05t°oM and n 



1 



+0.2 



These 



two independent constraints on R and Q are consistent. 

In summary we have attempted to assess the signifi- 
cance of the presence of acoustic oscillations in current 
CMB data using a model independent method. One 
could view Eq. ^ as a different class of phenomenological 
models which has the advantage (given the question we 
are attempting to answer) of smoothly interpolating be- 
tween the absence and presence of acoustic oscillations. 
We subsequently use the interpeak separation (or acous- 
tic oscillation frequency) to derive a constraint on the 
geometry of the universe. Such a constraint is indepen- 
dent to that derived from the position of the first peak. 
However it is necessarily included in any analysis which 
considers the full set of initial conditions put forward 
in 1^. With the rapid progress of CMB experiments, 
such an analysis will become essential to obtain accu- 
rate, model independent constraints on the cosmological 
parameters. 
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